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We consider information transmission through a noiseless quantum channel, where the information 
is encoded into massive indistinguishable partides: bosons or fermions. We study the situation in 
which the partides are noninteracting. The encoding input states obey a set of physically motivated 
constraints on the mean vàlues of the energy and particle number. In such a case, the determination 
of both classical and quantum capacity reduces to a constrained maximization of entropy. In the case 
of noninteracting bosons, signatures of Bose Einstein condensation can be observed in the behavior 
of the capacity. A major motivation for these considerations is to compare the information carrying 
capacities of channels that carry bosons with those that carry fermions. We show analytically that 
fermions generally provide higher channel capacity, i.e., they are better suited for transferring bits 
as well as qubits, in comparison to bosons. This holds for a large range of power law potentials, and 
for moderate to high temperatures. Numerical simulations seem to indicate that the result holds 
for all temperatures. Also, we consider the low temperature behavior for the three-dimensional box 
and harmònic trap, and again we show that the fermionic capacity is higher than the bosonic one 
for sufnciently low temperatures. 



I. INTRODUCTION 

The capacity of a transmission channel is an impor- 
tant quantity in several aspeets, ranging from fundamen- 
tal theoretical issues to technology; it is therefore inten- 
sively studied in classical [ÏJ and quantum Q information 
theory. The capacity of a channel depends on difïerent 
factors, like the type of information that one wants to 
transmit, the character of the physical material that re- 
alizes the transfer, the constraints that exist on the chan- 
nel, the assistance that is available for the information 
transfer, etc. 

The capacity of a channel for the transmission of classi- 
cal information by using quantum states is an important 
example, which has attracted a lot of interest recently. In 
this respect, a fundamental result, obtained more than 30 
years ago, is the "Holevo bound" @, i (see also @,[Í,0]). 
An essential message conveyed by the Holevo bound is 
that at most n bits (binary dígits ) of classical information 
can be carried by a quantum system of n distinguishable 
qubits (two-dimensional quantum systems). However, the 
Holevo bound includes arbitrary encoding and decoding 
strategies. This has the consequence that for information 
transmission with infinite dimensional systems through 
noiseless channels, e.g. with the modes of an electromag- 
nètic field, the Holevo bound prediets infinite capacities. 

A similar situation appears also for the capacity of 
transmitting quantum information, merely by definition. 

To avoid such nonphysical vàlues, it is important to 
introduce the relevant physical constraints, whilc max- 
imizing over the classical or quantum information that 
can be encoded in a physical system. Classical capacities 
of channels that carry photons in the form of modes of 



an electromagnètic field have been studied extensively in 
the last decade [!, d, [|[ . In this well as for the 

quantum capacity, as considered in this paper, the natu- 
ral choice is an energy constraint on the input states of 
the channel. 

In this paper, we will deal with channels carrying par- 
tides that have a non-vanishing mass. Such a problem 
is motivated by recent experiments producing atòmic 
waveguides in optical microstructures [ïo| . or in atom 
chips pH . These experiments indicate a new possibility 
of using massive partides in quantum channels for classi- 
cal as well as quantum communication over macroscopic, 
or at least mesoscopic distances. For such massive par- 
tides, it is natural to put an average particle number 
constraint along with an average energy constraint on 
the input state of the channel. We will show that with 
these constraints, the grand canonical ensemble of statis- 
tical mechanics attains the classical as well as quantum 
capacity of noiseless channels that carry trapped nonin- 
teracting massive partides. Moreover, we show that the 
classical and quantum capacities in the case of noiseless 
quantum channels (with the above physical constraints) 
are the same. Therefore in the remainder of this pa- 
per, the phrase "capacity of (quantum) channels carry- 
ing massive partides" (usually implicitly implying that 
the constraints on energy and particle numbers are al- 
ready taken into account), unless with the specific adjec- 
tive, will mean the classical as well the quantum capac- 
ity. We should mention here that sometimes, it is more 
natural to consider the average energy constraint with a 
fixed number of partides. In this case, the capacity is 
attained in the canonical ensemble. Calculations of dif- 
ferent properties of the canonical ensemble are usually 
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very hard, as compared to the grand canonical ensemble. 
However, several studies indicate that the two ensemblcs 
give similar vàlues for the average occupation numbers 
even for finite, but moderately large number of partides 
[Ï2| . However, there are physical quantities of interest 
that give drastically different vàlues. For example, the 
fluctuation in the average number of partides in a Bose 
Einstein condensate is unphysically large according to the 
grand canonical ensemble, while that in the canonical en- 
semble is physically meaningful 12]. In this paper, we 
will restrict ourselves to the grand canonical ensemblcs, 
for finding the capacitics of channcls carrying massive 
partides. The capacities, being dependent only on the 
average occupation numbers, will be similar already for 
a moderately large number of partides. 

Nonintcracting bosons exhibit a critical phenomenon, 
the Bose Einstein condensation (BEC). In this paper, we 
will show that the capacity of channels carrying nonin- 
teracting bosons also exhibit this criticality with temper- 
ature. The capacity changes its behavior from being con- 
cave to convex (with respect to temperature), at the crit- 
ical temperature. Noninteracting fermions do not show 
any critical behavior. However, interacting fermions can 
exhibit the Bardeen-Cooper-Schrieffer (BCS) transition. 
As we have shown in Ref. [Hj], the channel capacity in 
this case also indicates the BCS transition. It should 
be noted here that the convexity and concavity, that we 
speak about in this paper, are always with respect to 
temperature, and not with respect to mixing of states. 

An important motivation of this study is to compare 
the capacities of channels carrying bosons with that of 
fermions. We will obtain the following result analyti- 
cally: 

Capacities of noiseless quantum channels carrying non- 
interacting spinless fermions are higher than that of spin- 
less bosons for a wide range of power law potentials for 
sufficienly high temperatures. 

Since we obtain the above result by using perturbation 
theory up to the third order, it holds already for moder- 
ate temperatures. Note that our numerical simulations 
indicate that fermions are better carriers of information 
than bosons even in the low temerature region. 

It may be noted that although the quantum channel 
that we consider in the paper is noiseless, the states that 
are used to encode the information are allowed to be noisy 
(mixed states encoding). In the case of classical capacity, 
it will turn out that there exist ensembles of pure states, 
which can be used to encode the classical information (to 
be transferred) , for attaining the maximal information 
transfer. 

The paper is organized as follows. In Sec. QT] we discuss 
the Holcvo bound on classical information transmission. 
Wc definc the classical capacity of quantum channcls in 
the next section (Sec. IIIII) . In Sec. [TV] we define quan- 
tum capacity of quantum channels. In Sec. [Vi we discuss 
the important ensembles in statistical mechanics, and in 
its first three subsections, we discuss respectively the mi- 
crocanonical, the canonical, and the grand canonical en- 



sembles. In Sec. IV Dl we consider the situation when 
there is a constraint on the average number of partides, 
while there is no constraint on the (exact or average) 
number of partides. In the next section (Sec. IVI[) . we dis- 
cuss the capacities of quantum channels with constraincd 
inputs. We then briefly discuss the capacities of channels 
carrying photons in Sec. I VIII In Sec. IVIII1 we con- 
sider spinless noninteracting bosons and fermions, and 
calculate the capacity of channels in these cases. In Sec. 
IIX[ we prové a theorem stating that fermions can carry 
more information (classical and quantum) than bosons. 
In Sec. O we consider the low temperature case, for the 
3D box and the 3D harmònic trap. For both these cases, 
the fermionic capacity is again higher than the bosonic 
one for sufhciently low temperatures. In the last section 
(Sec. IXI[) . we make some concluding remarks. This pa- 
per presents on one hand the results of Ref. [l3j with 
more details, and on the other hand generalizes them to 
the case of trapped non- interacting bosonic and fermionic 
gases for a large class of power law potentials. Let us note 
that many of such potentials are feasible with currently 
available technology. We will see that the signature of 
BEC in the channel capacity in the case of a harmònic 
trap, is more pronounced than that in the case of a uni- 
form trap, as considered in Ref. [13j. 

It may be worthwhile to mention here that the mathe- 
matics required for obtaining the capacities, is similar to 
that in quantum statistical mechanics, where one maxi- 
mizes the entropy under different constraints. This is not 
very surprising, as the bounds on information transfer, 
discussed in Secs. [TT] and IIV1 are "entropy-like" quanti- 
ties. 



II. THE HOLEVO BOUND 

The Holevo bound is an upper bound on the amount of 
classical information that can be accessed from a quan- 
tum ensemble in which the information is encoded. Sup- 
pose therefore that a sender Alice (^4) obtains the clas- 
sical message i, and she knows that this happens with 
probability pt . She wants to send it to a receiver Bob (B) . 
To do so, Alice encodes the information i in a quantum 
state pi, and sends the quantum state to Bob. Bob re- 
ceives the ensemble {pi, pi}, and wants to obtain as much 
information as possible about i, for which he performs 
a measurement, that gives the result m, with probabil- 
ity q m . Let the corresponding post-measurement ensem- 
ble be {pi\ m , Pi\ m } ■ The classical information gathered 
can be quantified by the mutual information between the 
message index i and the measurement outeome pj: 

I(i : m) = H({ Pi }) - q m H({ Pllm }). (1) 

m 

Here H({r x }) = — J2 X r x l°g2 r x is the Shannon entropy 
of the probability distribution {r x }. Throughout the pa- 
per, we calculate the all the quantities on amounts of 
information transfer in bits (binary digits). Note that 
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the mutual information can be seen as the difference be- 
tween the initial disorder and the (average) final disorder. 
Bob will be interested to obtain the maximal information, 
which is the maximum of I(i : m) over all measurement 
strategies. This quantity is called the accessible informa- 
tion: 

hec = maxl(i : m), (2) 

where the maximization is over all measurement strate- 
gies. 

The maximization involved in the definition of accessi- 
ble information is usually hard to compute, and hence it 
is important to know bounds on I acc [H, 0, [Ü| • In partic- 
ular, in Ref. [1, 01, a universal upper bound, called the 
Holevo bound (or Holevo quantity), on I acc is given (see 
also 0, i ||): 

Iacc({pi, p t }) < x(ÍPi,Pi}) = S{p) - ^2píS(pí). (3) 

i 

Here ~p — ^íPipi is the average ensemble state, and 

5(ç) = -tr(çIog 2 ç) (4) 
is the von Neumann entropy of ç. 



III. CLASSICAL CAPACITY OF A QUANTUM 
CHANNEL: UNCONSTRAINED INPUTS 

Consider a quantum channcl 1Z that acts on d- 
dimensional quantum systems as inputs. Suppose that 
Alice wants to send some classical information i, that 
oceurs with probability pi, through this quantum chan- 
nel to Bob. She encodes this classical information in the 
quantum state pi, where the Hilbert space corresponding 
to the quantum states pi is d-dimensional. The classical 
capacity of this quantum channel is the maximal classi- 
cal information that can be sent through this channcl, 
and is therefore the accessible information of the ensem- 
ble {pi, IZ(pi)}, maximized over all such ensembles on the 
d-dimensional Hilbert space Ti d : 

C= max I acc ({pi,lZ(pi)}). (5) 

{pí,Pí} on H d 

However, capacities are usually defined in an asymptotic 
sense. Therefore, the classical capacity in this case is 

C°°=lim - max (6) 

n-*oa n { PuPi } on {H d )® n 



The Holevo bound implies that 



C < max 

{pi,Pi} on H d 



s(n(p))-j2piS(npi)) 



, (7) 



and 



lim — max 

n^oo n { Pi ,p,} on (H d ) i 



s(K(p))-J2 Pi s(n( Pi )) 



The quantity C°° is usually very difficult to handle. 
The Holevo- Schumacher-Westmoreland theorem [Ü], [Ï6[ 
states that in the particular case when the inputs are 
produets on the tensor product Hilbert space (Tl d )® n , 
the capacity (let us denote the capacity in this case by 

C T A ) is s iven b y 



max 
{pí-.Pí} on 7i d 



(9) 



For the case of the noiseless quantum channel that 
carries d dimensional quantum states noiselessly, all the 
above capacities equal to log 2 d. This value is attained 
by any complete orthogonal basis of pure states on TL d . 

Now for infinite dimensional quantum systems, the 
channcl capacity obtained in this way prediets an un- 
physical infinite capacity. This is because the Holevo 
bound itsclf docs not include any constraint on the avail- 
able physical resources in an actual implementation of 
the information transfer. In particular, arbitrary encod- 
ing and decoding schemes are allowed. To avoid this infi- 
nite capacity, one usually maximizes the accessible infor- 
mation over all ensembles that satisfies certain physical 
constraints. Due to the form of the Holevo bound on 
accessible information, such constrained maximizations 
are very similar to the ones in statistical mechanics. The 
same is true for the case of the quantum capacity which 
we briefly discuss in the succeeding section, and then in 
Sec. |V1 we briefly discuss some similar constrained max- 
imizations of statistical mechanics. 



IV. QUANTUM CAPACITY OF A QUANTUM 
CHANNEL: UNCONSTRAINED INPUTS 

We now consider the case of sending qubits (as opposed 
to bits) using quantum channels. The quantum capacity 
can be considered in (at least) the following four diffcr- 
ent situations [Ï7l ]: the quantum channel 7Z (acting on d 
dimensional quantum states at its input) 

• without the help of additional classical communi- 
cation (in this case, we call the quantum capacity 

Q°), 

• with an arbitrary amount of forward classical com- 
munication (Q - *), 

• with an arbitrary amount of backward classical 
communication (Q* - ), 

• with an arbitrary amount of both-way classical 
communication (Q* - *). 

Let us define the first case, the other definitions being 
similar. So, the quantum capacity Q° is 18] 



(8) 



Q° = sup lim 



log 2 D(n) 



(10) 
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where the supremum is over all such cases when there ex- 
ists a D(n) dimensional subspace §(n), of the total input 
space (7ï d )® n , satisfying the average fidelity criterion 

hm f (i/;\-R® n (\ï>){i}\)\i>)d\i>) = l, (11) 

where no classical communication was used in transfer- 
ring the input state from the sender to the receiver. 

In case of a noiseless channel, the fidelity criterion is 
automatically satisfied, and all the quantum capacities 
are equal to log 2 d. Now, if the dimension d of the input 
space is infinite, the capacities are again infinite, as in 
the case of classical capacity. 

There are several remarkable results that are known 
in the noisy case. In particular, Q° is the maximum 
coherent information [ÏJ, [ÏJ, Ül fèM '■ 

Q° = max \S(K(p)) - SÇL ® K($ p ))] , (12) 
p 

where the maximization is over all quantum states p de- 
fined on the Hilbert space TL , Q p is a purification of p, 
and I is the identity operator acting on quantum states 
of the ancillary Hilbert space that is required for the pu- 
rification. Furthermore [22l. [23|. 

Q° = Q^. (13) 

V. THE FUNDAMENTAL ENSEMBLES IN 
STATISTICAL MECHANICS 

In this section, we discuss the fundamental ensembles 
in statistical mechanics (see e.g. [2ÍL I25ll26j). 

A. The microcanonical ensemble 

Suppose that we have a physical system described by 
the Hamiltonian TL. We assume that the system has a 
fixed partició number N, and a fixed energy E. (To be 
more general, one must allow for fluctuations around E 
and N .) This represents a closed and isolated system. 
We want to find the state qmc of the system such that 
it maximizes the von Neumann entropy S. Let \E, N, k) 
dcnotc the state a energy E and particlc number N , and 
where k enumerates the degeneracy. Let íl(E, N) be the 
total number of orthogonal states with energy E and par- 
ticle number N, so that k — 1,2, ... , Çl(E, N). Since the 
system is isolated (as it has fixed energy and fixed number 
of partides), the set {\E, N, k)}^^'^ spans the allowed 
Hilbert space. Consequently, the maximum entropy is 
reached by the state 

Q(E,N) 

QMC = Sl(E,N) E \E,N,k)(E,N,k\, (14) 

which is actually the identity on the allowed Hilbert 
space. This is the microcanonical state of the system, 



and the ensemble {1/Q(E,N), \E,N,k)}^' N) is called 
the microcanonical ensemble. 



B. The canonical ensemble 

Consider next a physical system, described by the 
Hamiltonian TL, which has a fixed particle number N and 
a fixed average energy E. The average energy constraint 
forces every state g to follow 

ti(Hg) = E. (15) 

This physical system represents a closed, but not isolated 
system. We again want to find the state qc of the system 
such that it maximizes the von Neumann entropy S. One 
finds that [H 

É>c = ^exp(-/3W), (16) 
Zc 

the canonical state of the system. Here (3 = -j^r i with ks 
being the Boltzmann constant, and T the absolute tem- 
perature. The canonical partition function Zc is given 
by 

Z c = tr (exp(-/3W)) . (17) 

The state qc is called the canonical state of the sys- 
tem and the ensemble exp(— /3e,), is called the 
canonical ensemble, where {cí, \ei)} is the eigensystem of 
the Hamiltonian TL. Just like for the microcanonical en- 
semble, the particle number enters the calculations via 
the Hamiltonian and for determining the allowed Hilbert 
space. For example, in the calculation of the trace in 
Eq. (|17p . the summation runs over all combinations of 
different numbers of partides at different energy levels, 
under the constraint that the sum of all partides in all 
the levels is N. Moreover, for indistinguishable partides, 
we must also take care about the statistics of the parti- 
des. For example, if we have a trap, whose energy levels 
are e,, and in which N noninteracting spinlcss bosons are 
trapped, so that the bosons are described by the Hamil- 
tonian Tíb = X)i e i a l a i (where a\ and et^ are creation and 
destruction operators of the ith mode) , we have 

710=0,1,.. . ni=0,l s ... \ i ) 

(18) 

= E E ■■• II ex p(-^)' ( 19 ) 

no— 0,1,... ni— 0,1,... i 

v ■ ' 

E, nt=N 

where rii is the number of partides at energy level . For 
a given value of energy E, the temperature T is given by 
E=-JL{\og e Z c ). 
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C. The grand canonical ensemble 

The next step is to consider a physical system, de- 
scribed by the Hamiltonian Tí, which has a fixed average 
particle number N and a fixed average energy E. The 
average energy constraint forces every state g to follow 
Eq. (|15[) . while the average particle number constraint 
reads 



ix(Afg) = N, 



(20) 



where Aí represents the total particle number opera- 
tor. For example, for a system of spinless bosons, JV = 
X^i a í a íj where a\ and ai are the creation and annihi- 
lation operators of the zth mode. This physical system 
represents an open system. We again want to find the 
state qgc of the system such that it maximizes the von 
Neumann entropy S. One finds that the state is [Ü| 



Qgc 



■>gc 



■exp(-/3(«-/AA0)> 



(21) 



i.e., the grand canonical state of the system. Here fi is 
the chemical potential. The grand canonical partition 
function Zqc is given by 



Z a c = tx(exp{-0(H-nN))) 



(22) 



Note that in this case, the number of partides is not 
fixed, and in particular, 



oo oo 



^GC 



EE 

ng=0 n\=0 

11 i _ e -Kei 



exp {-^(ei-fS) < 



(23) 



in the case of noninteracting spinless bosons in a trap of 
energy levels {ei}, and where we have now denoted the 
chemical potential by /i fc . The ensemble 
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çrb 
^GC 



exp 



-/?E( £ « 



M fc )^J ,|ngX,...,n*,...)j 

(24) 

is the grand canonical ensemble, where the elements of 
the ensemble runs over all combinations of the n b, s {n\ = 
0,1,2,. ..,oo). 

As opposed to the case of the canonical partition func- 
tion, the grand canonical partition function requires the 
value of the chemical potential. We can determine the 
chemical potential from the constraint on the average 
particle number constraint, for given vàlues of the tem- 
perature and average particle number. We can then find 
the grand partition function with these vàlues. The av- 
erage occupation numbers, in the case of noninteracting 
spinless bosons, are given by 



1 



- 1 



(25) 



The average energy is given by E = J2i n i e i- Otherwise, 
for given vàlues of the average energy and average parti- 
cle number, the dual constraints on average energy and 
particle number can be simultaneously solved to obtain 
the chemical potential, and the temperature. 



D. The constraint on the number of partides 

One may now consider the case, when we do not have a 
constraint on the number of partides. In other words, ar- 
bitrary numbers of partides are allowed. For the energy, 
we still assume that the average energy is fixed to E. In 
this case, maximizing the entropy leads to the partition 
function 



oo oo 



Z» = £ £...exp 

71()=0 Tl±=0 \ i 

OO 

= n e exp (-P e i n i) > 

i m—0 



(26) 
(27) 



for noninteracting bosons in a trap of energy levels {et}. 
Note that this partition function is different from the 
canonical partition function in that there is no constraint 
on the number of partides (i.e. the constraint rii = N 
does not exist anymore), and from the grand canonical 
partition function in that the term P/J.AÍ is absent inside 
the exponent. 

In the case when the lowest energy level is zero or neg- 
ative, at least one of the sums in Eq. ([2T]) is divergent. 
Such divergence does not occur in Zc, due to the con- 
straint J2i n i — N. Neither does it occur in Zgc, as the 
term f3fiAÍ in the exponent "cures" the divergence. 

If the lowest energy level is positive, we have 



(28) 



and the average occupation number of the ith level is 

(29) 



" lo K ^= ' 



d(Pe t )- ae ~ e*«-l' 
so that the total average number of partides is 

Ar = E^ = E^ L· T' ( 30 ) 

i i 

while the total average energy is 

i 

Since all the energy levels are positive, these sums are 
convergent, and so we are not able to approach the "ther- 
modynamic limit" (N — > oo). The thermodynamic liniit, 
i.e. the limiting case of an infinite number of partides 
can be reached only in the trivial case when T — > oo. 
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This for example is the case for photons, in the form 
of modes of the electromagnètic field in a closed (finitc) 
cavity, in which case the zero frequency mode is not pop- 
ulated. With increasing cavity size, one may approach 
as close as possible to the zero frequency mode, but then 
the density matrix describing the system ceases, in the 
linht, to be trace-class. 

This problem does not arise in case of the canonical 
state, as in that case, if the lowest energy level is nonzero 
(positive or negative), we can define e° = e.; — e (eo is the 
lowest energy level) , so that we can rewrite the canonical 
partition function as 



Z b c = exp(-/3e iV) x 

E E 



exp 



"0=0,1, 



"1=0,1, 



Ei 

eyLp{-f3e N)Z, 



-N 



(32) 



where Z h Ca is the canonical partition function of a system 
of ./V noninteracting bosons in a trap with energy levels 
{e°}, and with the lowest energy level vanishing. The 
factor exp(— f3e N) in Eq. (f3"2")l is a constant for given 
particle number and energy. For the case of the grand 
canonical state, if the lowest energy is nonzero, we again 
rewrite the grand canonical partition function as 



^GC 



i "i=0 



— ■pr 



(33) 



where Zq Co is the grand canonical partition function of a 
system of N noninteracting bosons in a trap with energy 
levels {e^}, and with the lowest energy level vanishing, 
and with a correspondingly changed chemical potential 

& 

The discussion in this subsection so far concerned the 
case of noninteracting bosons. In the case of noninter- 
acting fermions, the Pauli exclusion principle forces the 
rii's in Eq. (|27[) to run from to 1, so that the parti- 
tion function in this case is (we consider "spinless" i.e. 
polarized fermions) 



of the value of the lowest energy level. Therefore, again 
we are not able to approach the thermodynamic limit, 
unless T — > oo. 



VI. CAPACITIES OF QUANTUM CHANNELS 
WITH CONSTRAINED INPUTS 

We now consider the case of quantum channels with 
constrained inputs. The motivation for such considera- 
tion is the fact that in the noiseless case, the capacities 
(both classical and quantum) are in some cases infmitc, 
and thus calls for putting some realistic constraints on 
the set-up. 

Surcly, the constraints will depend on the specific re- 
alization of the physical system under consideration. In 
anticipation of the cases that we will consider below, we 
put constraints on the average energy and the average 
number of partides of the input states. Suppose therefore 
that we have a system of partides that are sent through 
a channel 1Z, and that is governed by the Hamiltonian 
EL Let the corresponding total particle number operator 
be N. We only allow those inputs, p, to the channel that 
satisfies the following two constraints: 



and 



tr(HLo) = E, 



tr(Np) = N, 



(36) 



(37) 



for certain chosen vàlues of E and N. 

Consequently, the classical capacities C and C°°, are 
now respectively the same as the right-hand-sides of 
Eqs. © and ([ü]), but where the maximizations are over 
those ensembles whose average ensemble states satisfy 
the above constraints. For the case of C, the input en- 
sembles must satisfy exactly the constraints in Eqs. (|36|) 
and (f3"6"|) . For the case of C°°, the constraints must be 
suitably modified. In this case, they are (cf. Ref. (27j ) 



tr \ H í/°J = nE > 



(38) 



and 



Z * = II E exp (-P e i n i) 
= JJ(l + exp(-/3ei)). 



(34) 



The average occupation numbers in this case are (we are 
using the same notation for bosons and fermions) 



a 



log„ Z f 



1 



(35) 



so that the total number of partides and the total energy 
have fixed vàlues (for a given temperature), irrespective 



tr 



= nN, 



(39) 



where Hj = I ® I ® . . .W. . . ® I, i.e., it is the identity 
operator / at all the n sites, except at the ith one, where 
it is H. Similarly, N t = I <S> I <S> . . .N . . . ® I. 
It seems plausible that (cf. Eq. (O) 



TOC ,1 



(40) 



where the maximization is over all ensembles {pi, pi} such 
that the average ensemble state p satisfies Eqs. (|3"5)) and 
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(l37|) . The case of a noiseless channel with constrained 
inputs was considered in Ref. [27J ■ (Note that we are us- 
ing the same notations for the capacities as those in the 
unconstrained case.) In the succeeding sections, we will 
consider the right-hand-side of Eq. (|4D|) to be the "clas- 
sical capacity" of the noiseless quantum channel where 
the inputs are constrained by either or both the energy 
and number constraints. 

Similarly, it seems to be true that the quantum capac- 
ities (without, or with forward, classical communication) 
of a quantum channel with inputs constrained by Eqs. 
(0 and ©, are given by the Eqs. (fl"2")) and (flU)) respec- 
tively, but where the maximization in the right-hand-side 
of Eq. (fl"2"|) is over those inputs that satisfies Eqs. ([5]) 
and ((6]) . It is indeed true in the noiseless channel case, as 
we will show in Sec. IVIII Al (Actually, in the noiseless 
case, all the quantum capacities defined in Sec. IIVI are 
equal.) 



COMMUNICATION CHANNELS THAT 
CARRY PHOTONS 



VII. 



As we have already discussed in Secs. |HT] and IIV[ 
the classical and quantum capacities of a noiseless quan- 
tum channel that are constrained only by the dimension 
of the transferred state, predict infmite capacities. To 
avoid such infmite capacities, in the case of communica- 
tion channels that carry photons, one usually applies the 
energy constraint on the allowed states (or ensembles) 
that are fed into the quantum channel [E B II] • Suppose 
therefore that we have a system of photons that is being 
sent through a channel, and is described by the Hamil- 
tonian TL P = J2ifc UJ i a l a i- Here, we consider the case of 
classical capacity. From the considerations for massive 
bosons considered in Sec. IVIII Al it will be apparent that 
all the quantum capacties defined in Sec. IIVI are equal to 
the classical capacity, in this case of photons also. 

The energy constraint implies that for any ensemble 
{pi,Pi} used in a classical information transfer must sat- 
isfy 



tr (H p p) = E, 



(41) 



where E is a pre-assigned value of average energy for the 
system, and where ~p = ^ZíPíPí is the average ensem- 
ble state. The classical capacity in this case is therefore 
obtained by maximizing the Holevo quantity x({Pii Pi}) 
over all possible ensembles {pi,pí} that satisfy the aver- 
age energy constraint in Eq. (|4"T]) : 

(42) 



C P E = max XÍ{Pi,Pi})- 

{pí,Pí} satisfying Eq. (14 1 I t 



Therefore, 



Cl < 



{pi,Pi} satisfying Eq. I|411 



S(p). 



(43) 



However, the maximization on the rhs of the inequality 
(|4"3")l is the same as has already been considered in Sec. 



IV Dl since the number of photons is not conserved. The 
ensemble that maximizes the rhs of Eq. (|43|) consists 
of the number states (Fock states), which are mutually 
orthogonal. Therefore, the capacity C E is attained for 
this ensemble: 



c p 



H 



1 



exp 




(44) 



Here Z b is given by Eq. [[28]). and H({r t }) = 
— J^. ri log 2 Vi is the Shannon entropy of a set of proba- 
bilities {ri\. Rewriting the capacity in terms of average 
occupation numbers, we obtain the classical capacity of 
channels carrying photons as 



nP 



[rii log 2 m - (1 + Ui) log 2 (1 + m)] , (45) 



where n[ is given by the Eq. (|29|) . Although we have ob- 
tained this capacity by maximizing the Holevo quantity, 
the same result is also obtained by maximizing the ac- 
cessible information in the single-copy (nonasymptotic) 
case (see Eq. ([5])). 

The above classical capacity is obtained by using the 
average energy constraint. Several other types of con- 
straints on the energy are possible and all of them give 
the same classical capacity to lowest order (see Ref. [8j ) . 



VIII. CHANNELS THAT CARRY MASSIVE 
INDISTINGUISHABLE PARTICLES 

Let us now move on to the case of a communication 
channel that cames massive indistinguishable partides 
[Ï3j |. In this case, it is natural to impose the average 
particle number constraint along with the average en- 
ergy constraint. The classical capacity in this case is the 
maximum of the Holevo quantity over all ensembles that 
satisfy both thesc constraints, whilc the quantum ones 
are equal to the maximum of the entropy of the state 
that satifi.es both these constraints. As we will see, both 
classical and quantum capacities are equal for the noise- 
less channel. 

A system of noninteracting massive bosons exhibits 
a Bose Einstein condensation: Below a certain critical 
temperature, the number of partides in the lowest en- 
ergy level (ground state) becomes a significant fraction 
of what is present in the whole system. In Ref. [r|, 
we considered the channel capacity in Eq. (|5D|) below, 
and showed that it can exhibit the onset of the conden- 
sation. Here we present more details of the calculations 
presented there. The figures presented in Ref. [l3j were 
corresponding to the case where the bosons were confincd 
to a three-dimensional cúbic box with periòdic boundary 
conditions. Here we consider the harmònic trap, which 
is actually more close to usual experimental procedures. 
Unless mentioned otherwise, we consider in this paper a 
harmònic trap in three dimensions. 
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Let us note that noninteracting photons, despite be- 
ing bosons, do not show condensation, due to particle 
number nonconservation. However, effectively interact- 
ing photons may acquire an effective mass, and show 
condensation (see for instance Ref. (2^)- 



A. Noninteracting bosons 

Suppose that a quantum channel transmits spinless 
noninteracting bosons described by the Hamiltonian 
Tib = ^2i e t a l a ii an d that they are confined to a trap 
of energy levels {ei}. We will begin by considering the 
classical capacity. 



1. Classical capacity 

For a system whose average energy and average parti- 
cle number are given by E and N, the allowed ensembles 
{PíiPí} niust satisfy the average energy constraint 

tr(H b p) = E, (46) 
and the average particle number constraint 

tr(A/>) = N, (47) 
where J\í = aja^. The classical capacity given by 



Cbe 
E,N — 



Therefore, 



max 



{Pi,Pi} satisfying Eqs. (146 I t and 11471 



^E,N 



< 



max 



{Pí,Pí} satisfying Eqs. <46t and <47t 



x{{Pi,Pi})- 

(48) 

S(p). (49) 



Once again, the maximization of the rhs of inequality 
was considered before, in Sec. IV Cl There, it was 
discussed that the maximization of the rhs is attained by 
the grand canonical ensemble. The grand canonical en- 
semble consists of number states, which are mutually or- 
thogonal, and therefore, the channel capacity is attained 
for that ensemble, so that 



pbe tt 

^E.N — n 




exp -/?5>,- M b K 



(50) 



One may rewrite the capacity in terms of the average 
occupation numbers as (compare with Eq. (|45p ) 



Cbe 
E,N 



E 



log 2 1 



(51) 

Again, the classical capacity obtained here by maximiz- 
ing the Holevo quantity is the same as that obtained by 
maximizing the accessible information in the single-copy 
(nonasymptotic) case (see Eq. (|SJ)). 

2. Quantum capacity 



We now consider the case of quantum capacity. Con- 
sider the noiseless quantum channel that transmits ar- 
bitrary pure states of H® n , such that the energy and 
number of partides (and not their average vàlues) in any 
such pure state is nE and nN respectively. (Mixtures 
of such pure states are also noiselessly transmitted, by 
linearity.) For a noiseless channel, the average fidelity 
condition is automatically satisfied (see Eq. (fTTj) ). Con- 
sequently, the corresponding quantum capacity Q n , ir- 
respective of whether there is any classical side channel 
(see Sec. ITV)) . is given by 



Q n = log 2 n(nE, nN), 



where Çï(nE,nN) is the dimension of the subspace 
spanned by pure states of 7í®" having energy nE and 
number of partides nN. So, we are now in the micro- 
canonical ensemble. (To be more precise, one should al- 
low for small fluctuations around nE and nN, but the 
derivation is similar.) 

The system (described on the Hilbert space 7{® n ) con- 
sists of n identical parts, each being described on the 
Hilbert space 7í. Since the parts are identical, each part, 
on average, has energy E and number of partides N. 
Consider one such part, and let l(e,m) be the dimension 
of the subspace of that part whose energy is e, and which 
consists of m partides. (Again, we are disregarding fluc- 
tuations.) Then, 



log e n(nE, nN) = log e ^ l(e, m)to(nE - e, nN — ni ) 

e,m 

= log e J2 *(e> m ) ex P [ lo Se °(0 -l)E + E-e, (n -l)N + N- m)] , 

e,m 

i 



where we have 

E-t 
(n - 1)E 



< 1, 



N — m 
(n-l)N 



«1, 



as we suppose that n is large (compare with e.g. [29j|). 
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Therefore we can expand log e Í2((n — Í)E + E — e, (n — 1)N + N 
I 



m) as 



Therefore, 



log e fi((n - l)E + E - e, (n - 1)N + N - m) = log e íl((n - l)E, (n - 1)N) 



d log e ü 



dE 



(E-e) + dlog e ü 

E=(n-l)E,N=(n-l)N 8N 



E=(n-l)E,N=(n~l)N 



(N — m) + 



(52) 



\og e n(nE,nN) = log e íl((n -1)E, {n - 1)N) + log e ^ l(e, m)e 



Í3e+fi b /3m 



+ /3E- n b pN 



log e fi((n - l)E, (n - 1)N) + log e Z b GC (E, N) + f3E - n b /3N, 

I 



where Z GC (E, N) is the grand canonical partition func- 
tion of the system described on Tí, and having aver- 
age energy E and average particle number N. We will 
ultimately be interested in dividing the logarithm of 
Çl(nE, nN) by n and consider the limit as n — » oo, and 
in that case, the further terms on the right-hand-side of 
Eq. ([52"]) , after the third term, will not contribute. We 
now use the "fundamental equation for open systems" , 

-fc s Tlog e Z b GC = E - k B T{\og e 2)S(E, N) - (J'N, 

where S(E, N) is the von Neumann entropy of the grand 
canonical state of the system described on H, and having 
average energy E and average particle number N . (Note 
that we are using the information theoretic definition of 
entropy as in Eq. in which the multiplicative con- 
stant Ub is absent, and the logarithm is to the base 2.) 
Finally therefore, 

log 2 Q(nE, nN) = log 2 ü((n - 1)E, (n - 1)N) + S{E, N). 

Similarly, 

\og 2 n((n-l)E,(n-í)N) 

= log 2 Q((n - 2)E, (n - 2)N) + S{E, N), 

so that 

log 2 Q(nE, nN) = log 2 Q((n — 2)E, (n—2)N)+2S(E, N). 

This recursion can be carried on for p times, if n — p is 
large [30], so that we have 

log 2 Q(nE, nN) = log 2 Q((n-p)E, (n-p)N)+ P S(E, N). 
We moreover demand that lim — = 1, which can be 



satisfied even if n 
we have 

1 



p is large [31j . Since n — p <C n 



lim — log 2 — p)E, (n 

n^oo n 



■P)N) = 0, 



so that 



lim 



log 2 Vt{nE, nN) 



S(E,N). 



Therefore, the quantum capacity, Q b £ N , irrespective of 
whether we have an additional side channel for classical 
information transfer, of a noiseless quantum channel that 
transmits arbitrary states, p, of 7í, such that the average 
energy tr(W&p) = E and the average number of partides 
tx{Mp) = N, is given by 



Ïe,n= i™ 



\og 2 Çl(nE,nN) 



S(E,N). 



We therefore see that the quantum capacity is equal to 
the classical capacity obtained in Sec. IVIII A II 



3. Spinless bosons in a harmònic trap 

Let us now consider the case of noninteracting massive 
spinless bosons confined in a harmònic trap. The energy 
levels are therefore 



h(uj x n x + Lü y n y + Lü z n z ), n x ,n y , n z = 0, 1, 2, . 



(54) 



In evaluating the capacity numcrically, we must cut off 
the infmite sequence of energy levels at some point. The 
evaluated "capacity" of course depends on the point of 
the cut-off. In Fig. [TJ we plot the channel capacity for the 
case of 10 4 bosons (with respect to temperature) , where 
we show the "capacities" for the different points of the 
cut-off. In this case, the convergence is obtained roughly 
at about = 300. We see that the curves 

show a sort of fracture at a certain temperature, changing 
their behavior from being concave to being convex with 
respect to temperature. This can be seen more clearly in 
Fig. [2 where we plot the derivative of the plots in the 
preceding figure. For calculating the derivative, we use a 
four-point formula: 



dx 



1 



12A:r 

-f(x + 2/\x 



f{x ~ 2Ax) - 8f(x - Ax) + 8f(x + Ax) 
0(Ax 4 ). (55) 



In the thermodynamic limit, this bending causes the 
derivative to have a discontinuity, and the corresponding 
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FIG. 1: (Color online) Capacity (in bits or qubits) as a func- 
tion of the dimensionless variable T/T c for TV = 10 4 mas- 
sive bosons in a harmònic trap, with difïerent summation 
boundaries. T c denotes the critical temperature in the ther- 

.1/3 

modynamic limit. Therefore, T c = "" h ° ( y^rr I , where 



h^ha ( N \ 
k B \W) ) 

Lüho = (tüxUyLúz) 1 ^ is the geometrical average of the freqüèn- 
cies and Ç is the Riemann zeta function. For the purpose of 
all the figures in this paper, we consider a harmònic trap, and 
the case when The four curves (from bottom 

to top) represent cut-offs at = 40, 120, 220, 

and 300, respectively. Note that the last two curves almost 
coincide. 



temperature is the critical temperature of the condensa- 
tion. Fig. [3] shows the development of this discontinuity 
in the derivative of the capacity, with increasing average 
total particle number TV. The corresponding capacities 
are plotted in Fig. [U As seen in Figs. and UJ the 
critical temperature corresponding to the fracture grows 
with TV, as expected. The fractures indicate the onset 
of the condensation for the corresponding vàlues of N. 
The gap between the temperature corresponding to the 
fracture and the thermodynamic critical temperature is 
known to exist, and estimated vàlues of the gap has also 
been given (see Ref. [HI). Also, the capacities grow with 
the number of partides, as expected. 

It is known that the dimension of the system under 
consideration plays a role in determining whether a con- 
densation exists. For example, in the case of harmònic 
traps, the 3D and 2D traps exhibits condcnsations, while 
the 1D case does not show a condensation (see e.g. [33|). 
The cases considered in the previous (and latter) figures 
are 3D harmònic traps. In Fig. [ïïj we compare the quali- 
tative behavior of the derivatives of the capacities for 2D 
and 1D. 



FIG. 2: (Color online) Derivative of the capacity (in bits or 
qubits) as function of the dimensionless variable T/T c for TV = 
10 4 bosons. Note that the number of bosons is the same as 
that in the previous figure. Moreover, the curves in this figuro 
(from bottom to top) represent the cut-offs at 40, 120, 220, 
and 300, respectively, which are exactly the same as those in 
the previous figure. Note here also that the last two curves 
almost coincide. 



8000 




T/T 



FIG. 3: (Color online) Derivative of the capacity (in bits or 
qubits) as function of the dimensionless variable T/T c for dif- 
ferent number of bosons. Curves (from bottom to top) repre- 
sent 100, 500, 1000, 5000, and 10000 bosons, respectively. 
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FIG. 4: (Color online) Capacity (in bits or qubits) as func- 
tion of the dimensionless variable T/T c for difïerent number 
of bosons. At the critical temperature, the capacity changes 
its behavior from being concave to being convex (with respect 
to temperature). Again, the curves (from bottom to top) rep- 
resent 100, 500, 1000, 5000, and 10000 bosons, respectively. 



Going back to the case of the 3D harmònic trap, we 
consider energy as a function of the temperature, and 
see that the criticality is clearly visible in these curves. 
In Fig. [ül we plot the energy and its derivative with 
respect to temperature, for 500 spinless bosons in a har- 
mònic trap. Note that the capacity plotted with respect 
to energy does not indicate the condensation. 



B. Noninteracting fermions in a harmònic trap 

Let us now move over to the case of noninteracting 
fermions. Similar arguments as in the preceding subsec- 
tion imply that the capacities of N spin-s noninteracting 
fermions in a trap with energy levels {ei} are given by 

pfd _ n fd 



E,N 



H 



where 



2. 



■ exp 



f\ n f 

"i,mag 



GC 



mag=— s 



4c=n[ i+ex p(-^-^))] ! 



(56) 



(57) 



FIG. 5: (Color online) Derivatives of the capacities (in bits 
or qubits) are plotted as functions of the temperature for 
N = 100 bosons in one- and two-dimensional traps. The 
upper curve represents the 2D case, where the horizontal 
axis is of the dimensionless variable T/T^ D . And T% D = 



U(2) 



1/2 



where u>ï„ 



iíüxiüy) 1 ^ 2 33]. For the hgure, 
= uj y . In the one-dimensional 
case (lower curve), the horizontal axis is of the dimension- 

. It is 



less variable T/T* u , where T c 



we consider the case when uj x 



k B log c (2JV) i 

clear from the hgure that capacity of a one-dimensional trap 
does not show any onset of condensation, while that for a 
two-dimensional trap clearly does show it. 



magnètic quantum number. pJ represents the chemical 
potential in this case. The channel capacity in this case 
is reached by the fermionic grand canonical ensemble 



^GC \ i mag=- 



0i - ^ f )n[ 



mag 



f f f f f f \\ 

n 0,-si "0,-s+Ij ■ ■ ■ > n 0,si n l,~si ■ ■ ■ i n l,ai ■ ■ ■ I (i 



(58) 



where the elements of the ensemble runs over all com- 
binations of the n[ mag 's (n{ mag = 0,1, and mag = 
— s, — s + 1, . . . , s). Again the capacities can be rewrit- 
ten in terms of the average occupation numbers 



the power g = 2s + 1 being present due to the degener- 
acy of the spin states. The suffix "mag" represents the 



(59) 
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FIG. 6: (Color online) For N = 500 bosons, the dimensionless 
variable energy/ftw (where lj x = uj v = ui z — w) and its deriva- 
tive, are plotted against the dimensionless variable T/T c . En- 
ergy behaves quite similarly as the capacity, in that it changes 
its behavior from being convex to being concave, with respect 
to temperature, at the criticality. This is more clear from the 
derivative of the energy. 



FIG. 7: (Color online) For N = 100 fermions, the capacity (in 
bits or qubits) of spin-1/2 fermions is compared with that of 
spinless fermions. The horizontal axis is of the dimensionless 
variable T/Tf, where T/ is the Fermi temperature for spin- 
1/2 fermions. Note that we use the same Fermi temperature 
for both the curves. The Fermi temperature for a harmònic 

trap is given by T, = (t^- 



as 



níà _ 
E,N ~~ 



n fd _ 

— log 2 — 
9 9 




(60) 



In Fig. [71 we compare the capacities of spinless 
fermions with that of spin 1/2 fermions in a harmònic 
trap. 

Noninteracting fermions do not exhibit a condensation. 
However, interacting fermions can exhibit Cooper pair- 
ing, and consequently superfluid BCS transition. The 
channel capacity in such case also indicates the onset of 
the BCS transition (see Ref. [Ï3|). 



IX. FERMIONS ARE BETTER CARRIERS OF 
INFORMATION THAN BOSONS 

Numerical simulation with several vàlues of the total 
number of partides, with several types of traps, and over 
a large range of temperature reveals that fermions are 
better carriers of information than bosons. For sufh- 
ciently high temperatures, we have been able to obtain 
analytical results in this respect, as we have already pre- 
sented in the following theorem from Ref. fl3j ]. Here, we 



give a detailed proof of the theorem. We have already 
proven that the classical and quantum capacities are the 
same in the case under study. For definiteness, we will 
consider only the classical capacity in this section and 
the succeeding one. 

Theorem. For power law potential traps ( with power 7 
and dimension d ), and for sufficiently high temperatures, 
the capacity of spinless fermionic channel is better than 
that of spinless bosonic channel when 



1 1 1 

7 + 2 > 5· 



(61) 



Note that a power law potential with power 7 and 
dimension d is given by V (x±, x%, . . . , Xd) = r 1 , where 
r = y/ x\ + x\ + . . . + x 2 d , in the e?-dimensional Cartesian 
space {x\,x<z, . . . ,x d )- 

Proof. Let us start with the case of (spinless) bosons 
and perform the high temperature expansion. First we 
expand the fugacity 



z b = exp(^/i b ), 



in powers of N / S\ , where 



"exp(-/c/3ei), 



1,2,.. 



(62) 



(63) 
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Note that 

n = YA 

i 

= y 1 

i %b 

= y^Zb exp(-^e,) [1 - z b exp (-fie,)]" 1 
= z b S*i + z b 2 5 2 + z$S 3 + ■ ■ ■ ■ (64) 
Suppose now that 



N , ÍN\ 2 , /7V x3 



(65) 



Substituting z b , as in Eq. (|ü5|) , into Eq. and com- 

paring different powers of N / S\ , we find that 



bo = 0, h = 1, 

S*2 , 25f — 5i53 

& 2 = — ç-, &3 = ^ ■ 

Dl D-^ 



(66) 



We are now in a position to expand C b ^ N / log 2 e in powers 
of N/Si. We use the expression in Eq. ([S"Tj) for this 
purpose. In Eq. (|5ip. the substitutions in the different 
terms are done as follows. For the term n\ log e n b , for n b , 
we substitute 

n\ = z b exp(-/3e í ) [1 - z b exp(-/3e. i )] _1 (67) 
= z b e-^+z fc 2 e - 2 ^+z b 3 e - 3 ^ + .... (68) 



and for log 2 n b , we use Eq. ([ü7| to write 

log e n\ = log e z b - (3éí - log e (1 - z b exp (-/3ei)) 
= log e z b - /3éí 



(69) 



and this is then used for the substitution. For the term 
^1 + nfj log e ^1 + nfj , we use Eq. (|6"8"| to substitute for 

n b in both places, and then log e ( 1 + n\ ) is expanded in 



powers of n b . 

A similar calculation is done for spinless fermions. In 
this case, the fermionic fugacity 



Zf = exp (/3/r ) 

can be expanded as 

N ÍN\ 2 /iV x3 



(70) 



(71) 



in which the /i's are obtained from the average particle 
number conservation equation 

n = YM 

i 

= zf exp [1 + zy exp (-fie^)]" 1 

= 2/51-2^52 + ^53-... (72) 



as 



/o = 0, fx = 1, 

5 2 j, 25 2 — 5i53 



/2 — rT) /3 — 

Dl 



5? 



(73) 



Note that these are the same as in the case of bosons, 
except for the sign in / 2 . 

We perform the calculation up to the third order, and 
find that 



nbe 
^E.N 



E 



N 



.N . N 



a b A^ \ +&V log e — +fó [ — ) log 



log 2 e \SxJ 5i e 5i \Si J e 5i 



plus higher order terms, whereas 



A/V. . AT 
(74) 



C fd 3 / /v 



J 



N , N 



N 



lo g e 



AT 



plus higher order terms. The coefficients of first order 
perturbation are equal: 

al = a{ = Si + Dx, (76) 

where we have set 

D k = ^/fe I exp(-fc/?e I ), fc = l,2,.... (77) 

i 

In the next order, they differ by a sign: 



b _ / _ £2 



5 2 5 2 



" 2 = ^ 2 2 Si 
The third order perturbation coefficients are again equal: 

b / q o , "^3 _, 2 ^2 . 25| - 5i5 3 25 2 
a 3 = «3 = -^5 2 + — +— + -2 Dx — —D 2 +D 3 . 

o Dl Oi Dl 



--^Dx+D 2 . 



(78) 



Also, 



(79) 



ft = p{ = -5i, 



/3 2 b = 0, ^=25 2 . (80) 

Now, upto third order, the only coefficients that are 
different are those of {N/S^ 2 and (N/S^ 2 log e (AT/5i). 
We have 



= (3 b 2 < /3|. 



(81) 
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Also 



holds, when 



(82) 



We now evaluate ^ and ^ by using the density of states 



ç(e) = J d d pd d x S (j- + V(r) - e^j (83) 

for a system in d dimensionsal Cartesian space r = 
(x%, X2, ■ ■ ■ , Xd), trapped in a potential V(r). p — 
{P11P2, ■ ■ ■ iPd) denotes the respective momentum, with 
p — \p\. The integration in Eq. (|53")) is over the phase 
space. We may rewrite ç(e) as ïjji, 



ç(e) = constant • / dx\ . . . dxd{t — V) 

Let us now consider the potential as [3f 
V{r) = r 7 , 



(d-2)/2 



(84) 



(85) 



where r = \J x\ + . . . + x\. In this case, 



1 \ í + S-l 

ç(e) = constant- ei 2 



x / dïi . . . dïd ( 1 



so that 



/ \ d I d ~ 2 

ç(e) = constant • 2 . 



^2N7/2\ (d ~ 2)/2 



(86) 



where Xi = e~< (i = l,...,d). We are now ready to 
calculate the sums in Eq. ([52"]) . We have 



Dx J /3eexp(-/?e)ç(e)de 
~ /exp(-/3e)ç(e)de 

/ f3e exp(-/3e)eT +i ^de 

r / \ d I d ~ 2 

J exp(-/3e)e-< + 2 de 

r(í + ü=2) 

" 7 + 2- 



(87) 



The integrations over e are performed from to 00. Sim- 
ilarly, one can calculate and it is equal to + |). 
Therefore, 



Therefore, 



"2 < < a í 



(88) 



(89) 



1 1 1 



(90) 



That is, the fermionic capacity is greater than the bosonic 
one in such cases. However, we must now check for which 
potentials and dimensions, the above perturbation tech- 
nique is systematic. 

In the expansion of the capacity in terms of N/Si, 
given in Eq. (|74p for bosons, and in Eq. (|75p for fermions, 
the coefficients of £ , £ log e £ , (f) 2 , (f ) 2 log e £ , and 

(£) 3 , are all of the order /3~ ( -~ +2 - ) , since 



Si = exp(-/3e)e- ( T +i ^cZe 



/3" 



(91) 



and 



D i = J /3eexp(-/3e)e _( T + V)d e 



0" 



(92) 



Moreover, ~ p K ï ^' . For systematics of the expan- 
sion in Eq. (Í74|) for bosons, we need that the orders of (3 



a\f i: ftjj- log e £, c^(|-) 2 , and a|(f ) 3 should be in 
increasing order of p. This demand leads to the following 
condition: 



2rf , d d „ „ 
h<i>- + ->l>0. 

7 7 2 



(93) 



This requires that i + | > i, which is the same as the 
condition required for fermions having a higher capacity 
than bosons. Similar calculation for the fermions leads 
to the same requirement. 

Lastly, note that N/ S\ tends to zero implies that T — > 
00. This completes the proof. □ 

Remark 1. The condition in Eq. (fïïTj) includes e.g. the 
harmònic trap in 2D and 3D, the 3D rectangular box, 
and the 3D spherical box. 

Remark 2. In the proof, we work up to order (N/S1) 3 , 
and so the theorem holds for quite moderate tempera- 
tures. 

At this point, please note that we have numerically 
checked that the statement of the theorem holds also for 
low temperatures for the harmònic trap and for the 3D 
rectangular box [Ï3[ . In Fig. [5J we compare the capacities 
of channels carrying bosons and fermions for N = 100 for 
the case of a harmònic trap. 



X. LOW TEMPERATURE BEHAVIOR 

As we have already stressed, although the low tem- 
perature behavior of the capacity is not covered by the 
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FIG. 8: (Color online) Bosons vs. Fermions. We plot the 
capacities (in bits or qubits) of spinless bosons and spinless 
fermions with respect to the dimensionless variable T/T c , for 
100 partides. Note that the horizontal axis represents T/T c 
for both the curves. The lower curve is for bosons. The figurc 
clearly agrees with our theorem. Moreover, this figure along 
with other numerical simulations that we have performed, in- 
dicates that the theorem is true even for low temperatures. 
The inset compares the behavior in an exemplary case at high 
temperatures. 



Theorem in the preceding section, numerical simula- 
tions seems to indicate that the statement of the The- 
orem is indeed truc for lower temperatures. In this 
section, we want to point out that the statement can 
be proven analytically in some circumstances, for low 
temperatures. For example, for the case of a bounded 
volume V, containing N partides at temperature T, if 
the partides are bosons, then the capacity at low tem- 
peratures (specifically, for T < T c , where in this case, 

T c = h2 ^k fe) 1 ' is s iven by ( see e -§- ^ 

Cv,n,t = 5^(5/2) f mk B T \ i 
log 2 e 2^ 3 ^ 2vr ) ' 

i.c. the bosonic capacity ~ T* . On the other hand, if 
the partides are fermions, then the capacity, for T < Tf 

2 

(where in this case, T f = ( § W Í ) 3 )> is S iven b y 

(see e.g. [2o| ) 

^V,N,T TT 1 

log 2 e = TÏ>' 

i.e. the fermionic capacity seales as T. Clearly, the 
fermionic capacity is higher than the bosonic one for suf- 



ficiently low temperatures in a bounded volume, such as 
a 3D box. 

In the case of a harmònic trap in 3D, the bosonic ca- 
pacity seales as T 3 (see e.g. 37]), when the temperature 
is below critical. For sufficiently low temperatures, the 
fermionic capacity seales as T, which can be estimatcd 
by using the Sommerfeld expansions of Fermi functions. 
Therefore the fermionic capacity wins once again. 

XI. CONCLUSIONS 

In this paper, we have investigated the classical as well 
as the quantum capacity of noiseless quantum channels, 
carrying massive partides. We have considered spin- 
less noninteracting bosons and fermions. Noninteracting 
bosons exhibit Bosc Einstein condensation, and we have 
shown that this critical behavior can also be observed 
by studying the capacity of a quantum channel carrying 
bosons. We show that the capacity of such channels is 
concave (with respect to temperature) above the critical 
temperature, while it is convex (with respect to temper- 
ature) below that point. This criticality is absent in the 
case of fermions, as expected. We have numerically eval- 
uated the capacities of bosons and fermions for different 
numbers of partides. In the case of bosons, even in the 
case of a small number of partides, say 100, condensation 
can be observed from the qualitative change in behavior 
of the capacity. 

We have also shown analytically that a channel car- 
rying bosons is not as good a medium for transferring 
classical as well as quantum information, as a channel 
carrying fermions. This is true for a wide range of poten- 
tials that can be currently created in the laboratory. The 
analytical calculation for power law potentials holds for 
moderate temperatures. However, numerically we have 
checked that this is true even for low temperatures. It 
is tempting to believe that such superiority of fermions 
over bosons is genèric, at least for power-law potentials. 
In special cases, we have considered the low temperature 
behavior analytically, and have shown that for sufficiently 
low temperatures, the fermionic capacity is higher than 
the bosonic one, for the 3D box and the 3D harmònic 
trap. 
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